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Abstract 

We use a configuration interaction approach within the envelope function approximation 
to study the nature of the excitonic resonance in nano-hybrids, composite nanoparticles (NPs) 
combining a semiconducting and a metallic segment in contact. With reference to recent ex- 
perimental reports, we specifically study CdS-based nanorods with metallic NPs deposited at 
the tips (matchstick) or metallic coatings (core-shell). The excitonic states are computed taking 
into account both the renormalization of the electron-hole interaction and self-energy effects 
induced by the the metallic segment on the electron-hole pair, as well as by the dielectric envi- 
ronment, through an induced charge numerical approach. In neutral matchstick structures the 
metal NP has only a minor influence (~1 meV) on the excitonic states. When the metallic NP 
is charged the exciton becomes rapidly redshifted and spatially indirect. In contrast, in neutral 
core-shell structures the exciton energy redshifts by tens of meV. 

Introduction 

Hybrid nanoparticles -also called nano-hybrids- which combine material segments with different 
physical and chemical properties at the nanoscale open new venues to applications. For example, 
while quantum size effects allow to engineer very precisely the electronic properties of semicon- 
ductor nanocrystals with a number of possible applications, 1 combining semiconductor nanocrys- 
tals with a metallic nanoparticle 2 ' 3 - may result in a multi-functional nanomaterial, 4 with the metal- 
lic segment acting, e.g., as an electrical connection 5 or a preferred anchoring site for assembling 
of complex networks.-^ 

Nano-hybrids made from a semiconductor nanorod (SNR) with tips covered by a metal nanopar- 
ticle (MNP), 3 are one example of a growing set of nano-systems with optical properties which are 
not just the linear combination of those of the constituent segments. 8,9 A most important property 
of metal- semiconductor nano-hybrids, for example, is their potential for photo-catalysis. In CdSe- 
Au and CdS-Au samples the Fermi energy of Au is located in the semiconductor gap. After illumi- 
nation above the semiconductor gap, the electron-hole pair separates, and fast electron migration 



to the metallic sector takes place.— ^- As a result, the band gap recombination of the electron-hole 
pairs is almost suppressed, while the charged metallic tips are able to reduce a molecular acceptor 
in the solvent. Moreover, electrons can accumulate on the MNP until the Fermi energy aligns with 
the semiconductor gap, the number of accumulated charges being in the range of few units up to 
tens of electrons, depending on the dielectric constant of the solvent. 10-11 Such systems, based on 
the direct gap materials CdS and CdSe, absorbing in a wide part of the solar spectrum, are can- 
didates for efficient photo-catalysis, e.g., for hydrogen production. The mechanism for efficient 
charge migration has not been clarified, though. 

In an attempt to characterize the optical response of metal-semiconductor nano-hybrids, we 



shall discuss the nature of the excitonic states in systems with different geometries (see [figure][l][] 1 ), 
namely, a SNR with a MNP attached at one end, also called "matchstick" structure, and a SNR 
covered by a metallic shell (core-shell structure). We shall consider in particular situations in 
which the excitonic absorption taking place in the SNR segment is well separated with respect to 
the plasmonic absorption of the MNP or shell. This is the case of CdS-based structures, where the 
excitonic resonance can be clearly distinguished in absorption spectra.—*^ Interestingly, despite 
the exciton is uncoupled to the plasmon in CdS-Au samples, the exitonic absorption is still found 
to be blueshifted by few tens of nm in a matchstick structure with respect to bare CdS SNRs, in 
ensemble measurements. 12 
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Figure 1: Schematic cross section and 3D geometry of (a) a matchstick and (b) a core-shell metal- 
semiconductor nano-hybrid. 



Several energy terms contribute to the excitonic states in the nano-hybrids considered here. In 
a SNR the excitonic energy is primarily determined by radial quantum confinement effects. The 
longitudinal confinement, in turn, can be tailored continuously from OD to ID regimes. 14 In nano- 
hybrids, in addition to quantum confinement, carries also feel the self-polarization potential (image 
potential) induced at the semiconductor-metal interface. This is an attractive term for both electrons 
and holes in the region right outside the metallic surface. For example, in a MNP this may induce 
bound image potential states. 15 In a metal-semiconductor nano-hybrid, where a SNR and a metal 
are in contact, image potential and lateral confinement may compete, possibly affecting carrier 
localization. 16 Due to size quantization, the e — h interaction of photo-excited pairs is also expected 
to give a large contribution to the spectral properties of SNRs. Clearly, Coulomb interactions will 
be particularly effective in long SNR where the longitudinal confinement is weak. Furthermore, the 
overall interaction is strongly enhanced by the dielectric confinement induced by the low-dielectric 
constant environment, leading, e.g., to strong correlation effects. 17,18 The MNP induces a further 
renormalization of the e — h interaction in the opposite direction, tending to screen the interaction. 
Lastly, charging of the metal NPs after photoexcitation 1 1 will create an internal electric field which 
may strongly polarize the exciton pair, since it acts in the long, weakly confined direction of the 
SNR. 

In this paper we propose a modelization of excitonic states in metal-semiconductor nano- 
hybrids, focusing on CdS-metal systems with the matchstick and the core-shell geometries. Our 
modelization takes into account the electron-hole interaction exactly using the configuration-interaction 
method, with the semiconductor states within an envelope function. Electron and hole states at the 
excitonic resonance are assumed to be well confined in the SNR, 19 but the influence of the metal 
Coulomb screening on the exciton carriers is fully taken into account by a screening model. We 
will study, in particular, the relative contributions of the different energy terms discussed above, 
and the shift and nature of the excitonic resonance with respect to charging of the MNP. 



Method 

Our aim is to describe the excitonic complex in a SNR which is in close contact with a small MNP 
or a thin metallic shell, and surrounded by a dielectric environment. We assume that the plasmon 
resonance is well separated in energy from the excitonic one, and describe the influence of the 
metallic segment by a screening model. In other words, we calculate how the excitonic resonance 
is affected by the polarization of the metal due to an electron-hole pair excitation, neglecting the 
change in the polarization properties of the metal due to exciton population, which is reasonable 
in the linear regime. 8 The full Hamiltonian of the exciton in the nano-hybrid reads: 



ttNH _ tjSNR , zjpol . tjSNR , rjpol . ,Amp n , 

H X ~ H e,h + H e,h + H C + H C + V e,h ' (I) 



where 

iSNR ; 



Hyl is the single-particle Hamiltonian of the electron and hole in the bare SNR, neglecting 
the influence of the metal and the environment, 



tjSNR tjSNR , tjSNR ,r,s 

M e,h — H e + M h ■ \ z ) 



Here H s ff? is the single-particle Hamiltonian of the photo-excited electron (hole), which we 



e(h) 

describe with a 3D single-band envelope function formalism, 

with i = e,h. p is the momentum operator, m* the effective mass and V c the step-like con- 
finement potential due to the different band alignment of the semiconductor with the metal 
and the insulating environment. 

H p e ° h represents the influence of the metal and the environment on the SNR single-particle 
states, that is, the self -polarization (image) potentials: 



H$ = V s P(r e )+V s P(r h ). (4) 

Details about the calculation of self-polarization and polarization terms are given in Refs.— "21 
See also Supporting Information. 

Hg R is the usual electron-hole Coulomb interaction term for homogeneous dielectrics, in- 
cluding local screening: 



Hi? R (T e ,T h ) = - T ^— , (5) 

c-S I *e *-h I 

where £ s stands for the semiconductor dielectric constant. 

Hg° is the Coulomb attraction between the hole and the electron polarization charges in- 
duced near the interfaces. The term, which will be projected onto a basis of single-particle 
exciton states, can be written as: 



fl?V«,r*)=£|r)4^H (6) 

r,s \ r e r h\ 

where \r) and \s) are single-electron eigenstates, while p l rs (r e ) is the charge density induced 
by the electron density charge \r)\s) at the position r e . See Supporting Information for further 
details. 

V™ p is the single-particle potential exerted by N electron charges trapped in the MNP on 
the photo-excited electrons and holes. We model trapped electrons as a multiply-charged 
impurity located at the center of the MNP. Then, V^ p = Vi mp (r e ) — Vi m p{^h), with: 



Vim P (r) = -^—-+V^(r). (7) 

Here tq is the position of the impurity and £ m the dielectric constant of the metal. V?° is 



the Coulomb interaction with the polarization charges induced by the impurity on the metal 
surface, 

rP° l („\ __ f PN\ r 0> 



J l r o r l 



In an ideal metal e m = °°, so the first term in Eq. ([equation] [7] []7) is zero and all the in- 
teraction takes place through the surface charges. However, owing to numerical limitations, 
in our calculations we shall consider a quasi-metal with large but finite e m (see Supporting 
Information). Yet, the dominant contribution is by far that of Vf°, so that we capture the 
essential features of the MNP screening. 

Note that we have assumed the exciton feels static screening from the metal. This is an accu- 
rate description for nano-hybrids where the excitonic resonance does not overlap with interband 
transitions of the metal, which is the case for most metals and many nano-hybrids of interest.- 12 ^ 3 - 
For CdS-Au nano-hybrids, however, gold interband transitions are already present at the frequency 
of the excitonic resonance, 12 which has an impact on the dielectric response of the metal.— In 
this case, a static screening model implies assuming that the metal is relaxed when we probe the 
exciton. This is likely to be the case for photoluminescence spectroscopy, but should be taken with 
caution when comparing with absorption measurements. Also, for CdS-Au the spectral overlap 
between semiconductor and metal excitations implies Forster energy transfer may be important 
in describing the exciton dynamics. 25 Indeed, this mechanism has been proposed to explain the 
photoluminescence quenching observed in CdSe-Au nano-hybrids. 3 In CdS-Au systems the influ- 
ence is less clear because Forster processes are weaker when excitons and surface plasmons are 
decoupled.— 

Starting from a cylindrical SNR with radius R and length L, below we shall study: i) matchstick 
structures ( [figure][l][]ir a)) formed by attaching at one end of a SNR a metal hemisphere of radius 



R m , and ii) core-shell structure ([figure][l][]l [b)) formed by coating the semiconductor by a metal 



shell of thickness AR. Both the SNR and the metallic sectors have rounded corners of radius Rq. 



A finite curvature is needed to avoid numerical instabilities and is consistent with TEM images.— 
Note that the exact shape of the MNP is not known. Our results do not depend sensitively on the 
exact choice of the MNP geometry, unless otherwise stated. 

The nano-hybrids are supposed to be surrounded by an insulating environment extending to 
infinity. Typical values for the confinement potential exerted by the dielectric environment are in 
the few eV range. The metallic layer has a Fermi energy which lies in the gap of the semiconductor 
and has a continuum of states at the semiconductor conduction edge. Nevertheless, in general 
the semiconductor states at the gap do not hybridize easily with the metallic continuum, and are 
well confined in the semiconductor. 26 For the specific Au/CdS interface it has been shown in 
Ref. 19 that the energy bands near the gap are very little affected by Au adsorption on CdS. In 
order to incorporate this effect in the envelope function approximation at the single band level, 
we include in our model a large barrier for the semiconductor states at the semiconductor-metal 
interface. For simplicity, we use the same barrier height as for the semiconductor-environment 
interface. Although the barrier height and precise position are to some extent arbitrary, we have 
checked that moderate variations do not affect our results qualitatively. To model the carriers in 
the semiconductor, we use CdS material parameters (see Supporting Information). 

Results and discussion 
Matchstick structure 

CdS nanorods can be grown with ultra-narrow diameter and relatively narrow distribution, from a 
few units to below the Bohr radius,— 12 -^ and high aspect ratio (from tens to hundreds nm). Here, 
we take R — 2.5 nm and selected lengths in the range L = 5 — 20 nm. In this case the diameter 
slightly exceeds the Bohr radius. We assume the MNP radius is R m = 2.5 nm. 



Neutral MNP 

We start by investigating the different energy contributions with a neutral MNP. In this case, V^ p 
and the Hamiltonian reads: 



jNH 



Hx =H e ,h + Hc, 



(9) 



where H e i l = H^ R -\-H^° h are the single-particle terms and He = H^ NR +Hfj° the electron-hole 
interaction. The complex dielectric configuration influences the excitonic resonance via the image 
potential H v e ° h and the Coulomb polarization term, H^° . To understand the role of these terms, in 
[figure][2][]2f a) we show the calculated image potential along the longitudinal axis of the struc- 
ture. Carriers are confined inside the SNR, whose dielectric constant is intermediate between that 
of the dielectric environment and that of the metallic tip. Therefore, the potential V sp felt by photo- 
excited electrons and holes has (i) a repulsive character (V sp > 0) in most of the semiconductor 
region, owing to the influence of the dielectric environment in the radial direction, (ii) a poten- 
tial barrier at the interface with the dielectric environment (z ~ 4 nm), and (iii) a potential well 
near the interface with the MNP (z ~ 16 nm). A more complete perspective of V sp is given in 



[figure][3][]3f a), which shows the corresponding 2D plot. 
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Figure 2: Self-polarization potential along the longitudinal axis of (a) a matchstick structure with 
L = 12.5 nm and (b) a core- shell structure with L = 7.5 nm. The dashed line highlights the zero 
potential position. Background colors indicate the type of material in each region: semiconductor 
(brown), metal (yellow) or dielectric (green). 

Next we consider the Coulomb polarization terms, H^° . For a preliminary insight into its influ- 



ence, [[figure][3]n3fb) shows the polarization charges induced by the electron in the lowest-energy 

9 



orbital (\r) = \s) = |0)), p' e0 , which is the dominant term in Eq. ( |[equation][6][]6[ ). One can see that 
the charges distribute near the interfaces, and they have different sign on the semiconductor inter- 
faces with the environment (where they are negative) and metal (where they are positive). Thus, 
Coulomb interaction with the hole will be enhanced (weakened) near the dielectric environment 
(conducting metal). 



O 






O 






o 
U 




Figure 3: Contour plots of different terms arising from dielectric mismatch in a matchstick (upper 
row) and a core- shell structure (lower row), (a) and (d): self-polarization potential, V sp . (b) and (e): 
polarization charge induced by the dominant electron configuration, p l e0 . (c): Coulomb potential 
arising from N = 2 electrons localized in the center of the MNP, Vi mp . 

We are now in the position to perform quantitative estimates of the exciton energy in matchstick 
[figure] [4] []4[ a) shows the single-particle contribution to the energies of an electron- 



structures 



hole pair as a function of the SNR length. To highlight the different contributions, we compare 
a bare SNR system (no dielectric mismatch effects, red squares) with a SNR in a dielectric envi- 
ronment (blue circles) and with the attached MNP (green triangles). In all cases the confinement 
energy decreases with increasing L due to the relaxation of the longitudinal confinement.—^ In 
addition, the dielectric environment is responsible for a large, rigid blueshift of about 100 meV, 
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independent of L. This is due to the repulsive nature of V sp , mostly arising from the radial di- 
electric confinement induced by the insulating environment. The inclusion of the metal has the 
opposite effect, namely to redshift the electron-hole pair energy, due to the attractive potential near 
the metal interface. However, this contribution is only significant for short L. For structures with 
typical length of L ~ 15 nm, it is about 1 meV. 

A similar trend holds for the exciton energy, i.e., including Coulomb interaction terms, as 



shown in [figure] [4] []4^b). There are however two visible differences. First, the magnitude of 



the blueshift due to the dielectric environment is now reduced to less than 20 meV. This is because 
H^° shifts the energy in opposite direction with respect to V sp ,2&Q with the the latter being slightly 
larger than the former, 18 which leads to a partial cancellation. Second, the blueshift decreases for 
short SNR. This is because Coulomb interaction terms, including H^° , are enhanced in strongly 
confined structures, while V sp remains as a short-ranged interaction near the interfaces. 
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Figure 4: Energy contributions to the excitonic energy of a matchstick structure, (a) single-particle 
energy of the electron-hole pair, neglecting dielectric confinement effects (red squares), includ- 
ing screening from the environment (blue circles) and that from environment plus metal (green 
triangles), (b) same but for excitons, including Coulomb interaction terms. 



One important point concerns the nature of the excitonic transition. In particular, one may 
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wonder whether the excitonic wavefunction is distorted or even localized by the presence of the 
metallic segment and the ensuing image potential well. To this end, we calculate the average 
electron and hole charge density, taking the charge density of the exciton ground state and inte- 
grating over the other particle's degrees of freedom. For example, the electron charge density in 
the excitonic ground state is evaluated as 

Pe(r e )= [Vx(T ei T h )*V X {Te,T h )dT h (10) 



where ^(r^r/,) is the lowest-energy eigenstate of Eq. ( [equation] [1][]1 ). An example is shown 



in 



[figure] [6] [] 6 (N = 0) for a structure of intermediate length. Both the electron and hole are 
localized in the middle of the SNR, as would be the case in the absence of the MNP. No sizable 
distortion is observed, despite the asymmetry of the structure. As we shall also discuss for core- 
shell structures, the image potential is usually too short ranged to localize carriers in CdS. However, 
this also depends on the effective mass and the quantum confinement. For specific materials, 
such as GaP, 21 and/or crystallographic direction, therefore, one might obtain interface localized 
excitons. 

Charged MNP 

Next we investigate the influence of charging the MNP, which might be the result of fast electron 
migration to the metallic sector and hole scavenging.—^ Upon continuous optical excitation, a 
new exciton is formed in the presence of the electron(s) trapped in the MNP. The Hamiltonian is 
then given by: 

Hg H = H eth + V™ p +H c . (11) 

The electrons trapped inside the MNP redistribute over the metal surface yielding a potential 
Vi m p, as plotted in [figure] [3] [] 3 c). The effect of this additional term is shown in [figure][5][]5[ 



with the MNP tip charged with up to three electrons. The figure shows the (a) electron confinement 
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energy, E e ; (b) hole confinement energy, £/,; (c) electron-hole binding energy, E^, and (d) exciton 
energy E x =E e + E h +E b . 
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Figure 5: Influence of N electrons trapped in the MNP on the energy of carriers confined in the 
SNR. (a) electron energy, (b) hole energy, (c) exciton binding energy, (d) exciton energy. 

Being V ;mp a repulsive term for electrons, it shifts upwards the electron energy, E e , as shown in 



[figure] [5] [] 5 'a). The repulsion is weaker for the larger L because the electron localizes on the far 
side of the SNR to minimize the repulsions with the metal. As a result, E e is more sensitive to L in 
the presence of a charged MNP than with neutral MNP (compare, e.g., N = 3 and N = curves in 
panel (a)). 



For holes, see |[figure][5][]5| ^b), Vj mp is attractive. Accordingly, the confinement energy, E h , 
becomes negative upon charging, and the L dependence almost disappears because holes become 
trapped at the semiconductor-metal interface, which clearly dominates over the weak longitudinal 
confinement. 



[figure][5][]5"c) shows the electron-hole binding energy, Ej,. Its absolute value decreases with 



N, since electron-hole pair is pulled apart by the electric field induced by the charged MNP, reduc- 
ing their mutual Coulomb interaction. An interesting point is that this effect tends to saturate at 
N = 2, since at this value the electron-hole pair is already as splitted as the longitudinal confine- 
ment allows. An additional charge (N = 3) does not substantially change the binding energy for 
the long structures. The Af-dependence remains efficient only for the shorter structures, where the 
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quantum confinement forces the electron and hole to share the same spatial region. 



The exciton energy, E x , see |[figure] [5] []5f d), turns out to be dominated by the hole attraction 
due to the charged MNP. As a result, the exciton energy decreases with N, in a similar fashion to 
holes. Yet, excitons are more sensitive than holes to the rod length, which is a consequence of the 
significant variations of E e and E b with L. For long structures, the shift is a large fraction of an eV, 
and the resonance may fall in the broad plasmon resonance. 12 

The charge densities as a function of the number of electrons in the MNP are plotted in 



[figure] [6] [] 6 For each N, left (right) panels show the electron (hole) charge density. It can be 
clearly seen that an increasing electron charging pushes electrons and holes in opposite directions, 
despite the reciprocal Coulomb attraction. A striking feature is that, for N > 2, holes do not simply 
localize over the entire semiconductor/metal interface. Instead, they show localization in a ring- 



like state. This result can be understood from the shape ofVj mp . As can be seen in [figure] [3] [] 3 c), 
Vimp is stronger on the metal/dielectric interface than on the metal/semiconductor one, because of 
the larger dielectric mismatch. Besides, it is stronger near the corners of the hemisphere owing to 
the larger curvature, which translates into a larger gradient of dielectric constants. This favors hole 
localization near the corners of the SNR. 32 Note that holes are particularly sensitive to dielectric 
mismatch effects (an attractive potential in this case) in the radial direction because of the larger 
in-plane mass, m^(_L), as compared to that along the growth direction m* h (z). 

Core-shell nanorods 

In this section we study the effect of the dielectric mismatch in a SNR with R = 2 nm, L = 7.5 nm, 
coated with a metal shell of thickness AR = 0.5 nm and embedded in a dielectric environment. 



[figure] [2] [] 2^b) shows the longitudinal profile of the image potential. One can immediately 
spot important differences as compared to the case of the matchstick structures, namely, (i) V sp is 
now attractive in the core, and (ii) there are only potential wells near the semiconductor interfaces. 
The differences are due to the dominant contribution of the metal shell all over the SNR, as op- 
posed to the matchstick case, where most regions felt a dominant contribution from the dielectric 
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Figure 6: Localization of the electron (e) and hole (h) of an excitonic pair in a matchstick with 
L = 12.5 nm as a function of the number of electrons trapped in the MNP. 



environment. The radial profile of V sp is similar to the longitudinal one, as can be seen in the 2D 



plot of [figure] [3] [] 3 d). The polarization charges induced by the electron in the lowest-energy 
orbital, p l e0 , are illustrated in [figure][3][]3| ^e). By comparing with the results for corresponding 
matchstick structure, panel (b), one can also notice qualitative differences. Namely, the charges 
are now repulsive for the hole all over the semiconductor surface. 

The effect of the polarization terms on the energy levels is studied in [figure][7][]7f a). Solid 
(dashed) lines are the results including (excluding) such terms. For the single-particle electron- 
hole pair (thin lines), the image potential translates into a redshift of AE e ^ = —56.7 meV. For 
the interacting exciton (thick lines), this is partly quenched by the polarization charges weakening 
Coulomb interaction, so that the redshift reduces to AEx = —41.5 meV. 



V sp in |[figure][3][]3f d) suggests that carriers in core-shell structures may localize in the poten- 
tial well near the metal surface. The large radial mass of holes in CdS, m*(_L) = 2.56mo, would 
further support this phenomenon because it entails moderate kinetic energy in spite of the strong 
radial confinement. The situation may however change for excitons, because H^° is expected to 
counteract V sp . To investigate this point, we have calculated the charge densities of the electron 
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Figure 7: Influence of the dielectric confinement in a core-shell nano-hybrid. (a) single-particle 
energy of the electron-hole pair (thin lines) and excitonic energy (thick lines), (b) and (c): radial 
cross-section of the electron and hole charge density in the excitonic complex, respectively. Solid 
and dashed lines include (exclude) dielectric mismatch terms. In (b) and (c), dotted lines are used 
to represent single-particle densities. 



and hole within the exciton complex and plotted them in [figure] [7] []7 b) and (c), respectively. 
Solid (dashed) lines include (exclude) polarization terms. One can see that the inclusion of such 
terms visibly reduces the density in the center of the rod, but no drastic localization near the metal 
occurs. As a matter of fact, the single-particle states (dotted lines) do not localize either. This 
shows that the image potential well formed by V sp is generally too thin to compete against quan- 
tum confinement. 



Conclusion 

We have developed a model to study excitonic resonances in semiconductor-metal nano-hybrids, 
whereby the excitonic states are subject to the screening of nearby metallic segments and dielectric 
environments. We have explicitly estimated the influence of the different energy contributions to 
the excitonic frequency in the case of CdS -based nano-hybrid with the matchstick and the core- 
shell geometry. 
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In neutral matchstick structures, the metal sector has a minor influence on the energy and 
charge density localization. This suggests that the blueshift observed in related experiments upon 
Au absorption 12 does not arise from polarization effects, rather from structural changes or by in- 
teractions with the interband transitions of the metal. The influence of the metal becomes however 
critical when the MNP is charged. The exciton becomes rapidly redshifted with increasing number 
of charges in the metal. Moreover, the Coulomb interaction with the additional charges overcomes 
the electron-hole interaction, and the exciton becomes spatially indirect. The sensitivity of the ex- 
citonic resonance upon charging of the MNP suggests the use of matchsticks for sensing of charged 
molecules. 

In core- shell nano-hybrid, on the contrary, the uniform metal coating redshifts the exciton 
energy by tens of meV. The stronger influence of the metal as compared to the matchstick is due 
to dominant role of the radial confinement. For certain nano-hybrid geometries, we predict the 
existence of exciton states where the hole is attached to the lateral sides of the SNR. 

Associated Content 

Supporting Information. Details on the method of induced charges computation. Calculation 



procedure and parameters used for solving Eq. ( [equation] [1] [] 1 ). This material is available free 



of charge via the internet at http://pubs.acs.org 
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Supporting Information 

Induced charge computation method 

Here we summarize the theory of the screening model we use. We start from the Poisson equation 
for inhomogeneous dielectric media, 

V[e(r)Vvr(r)] = -4wp(r), (12) 

where e(r) is the position-dependent dielectric constant, \j/(r) the electrostatic potential and p(r) 



the source charge (electron, impurity, etc). We want to replace Eq. ( [[equation] [1 2] [] 12] ) by the 
Poisson equation in the vacuum of the source charge plus an additional induced charge, p;(r), 
which implicitly accounts for local screening and interface effects, 

VV(r) = -4jr(p(r) + p;(r)). (13) 

Taking Eq. ( |[equation][13][]13| ) into Eq. <\ [equation] [1 2] [] 12] ) one obtains: 



The first term of p,(r) is the contact interaction correction, pf. When added to the source charge, it 
gives the locally screened charge, p(r) + p ; c (r) = p(r)/e(r). The second term is finite only when 
Vg(r) 7^ 0, and it gives rise to the interface effects. We can then split the electrostatic potential in 
Eq. ( |[equation][13][]13[ ) as \j/(r) = y^(r) + W P (r), so that: 



7l„,-S ,-PW 



V 2 ^(r) = -47T^, (15) 

e(r) 



72 



V> p (r) = -Anpi. (16) 



where p = p, — pf is the induced charge excluding the contact charge (first term of [equation] [14] [] 14 1. 
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tyAr is the usual homogeneous electrostatic potential, as in Eq. ( [[equation] [5] []5[ ), while \j/ p is the 



polarization potential, as in Eq. ( |[equation][6][]6[ ). Further details about the calculation of p,- can 
be found in Ref. 20 



Calculation procedure 



Hamiltonian ( [equation] [1][]1 ) is solved following the scheme of [figure] [8] [] 8 



Image potential 




le-particle states (e,h) 



ngte-p, 



Coulomb integrals 



\n 



Klomb polarization 
integrals 



Figure 8: Flow diagram showing the procedure for solving Eq. ( |[equation][l][]l| >. 

First, knowing the geometry and dielectric constants of the system, we calculate the image 
potential, V sp . This is done numerically using the induced charge computation method for arbi- 
trary geometries with cylindrical symmetry. 20 Solving Eq. ( [equation] [ 1 4] [] 1 4| ) for the self-energy 
requires a 3D discretization of the space, with a high density of points near the metal interface, 
where induced charges are essentially localized on the surface. This leads to huge and dense linear 
systems Ax = b, which we solve using a recently developed library combining LU decomposition 
of the matrix A with state-of-the-art out-of-core algorithms. 21 If the MNP is charged, we further 

rPOl 



compute the potential induced by the electrons near the metal surface, V t , Eq. ( |[equation][8][]8[ ) 



The induced charge p' N (r' ) is also calculated using the induced charge computation method. 

Second, we determine the single-particle electron and hole states. To this end, we integrate 
Hfi R +Heh +^T/j P usm g cylindrical coordinates and finite differences. We then build a basis 
set to expand the exciton states as Hartree products between electron and hole states. The basis 
includes 20 low-energy electron spin-orbitals (12 with azimuthal angular momentum m z = and 4 
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with m z = ±1), and 40 low-energy hole spin-orbitals (24 with m z = and 8 with m z = ±1). Next, 
the Coulomb interaction is taken into account, H^ NR + H^,° . p l rs is obtained following Refs.,— ^- 
and the Coulomb matrix elements are evaluated using Monte-Carlo routines. Finally, the Hamilto- 
nian matrix including single-particle states and Coulomb matrix elements is diagonalized using a 
configuration interaction method,— - — which yields the exciton states. 



Parameters 

The same material parameters are used for both matchstick and core-shell structures. We use 
CdS effective masses. For electrons we take m* = 0.21 mq, where mq is the free electron mass. 36 
The hole effective mass in wurtzite structures is strongly anisotropic. We take mt(z) = 0.22niQ 
(m^(_L) = 2.56mo) for the longitudinal (radial) direction. These values correspond to the mass of 
C-band holes, which is the dominant component of low-energy holes in CdS nanorods. 37 V c is zero 
inside the semiconductor and 4 eV outside (the electroaffinity of the semiconductor). The dielectric 
constant is e s = 9.2 in CdS, 3 £ e = 3.0 in the dielectric environment and S m = 100.0 in the metal. 
Increasing further the dielectric constant of the metal has only a minor quantitative effect on the 
results, at the price of a fast increase of the computational burden. A rapid but continuous variation 
of e in a region of 5 A around the interfaces is assumed to avoid singularities in the self -polarization 
potential. 38 In all cases we take a small corner radius, Rq = 0.25 nm (see [figure] [!][]![ )■ 
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